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Abstract 

We introduce higher order polynomial deformations of A\ Lie algebra. We con- 
struct their unitary representations and the corresponding single- variable differential 
operator realizations. We then use the results to obtain exact (Bethe ansatz) solu- 
tions to a class of 2- mode boson systems, including the Boson-Einstein Condensate 
models as special cases. Up to an overall factor, the eigenfunctions of the 2-mode 
boson systems are given by polynomials whose roots are solutions of the associated 
Bethe ansatz equations. The corresponding eigenvalues are expressed in terms of 
these roots. We also establish the spectral equivalence between the BEC models 
and certain quasi-exactly solvable Schor dinger potentials. 
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1 Introduction 

Polynomial algebras are non-linear deformations of Lie algebras and have recently found 
widespread applications in theoretical physics whereby they appear in diverse topics such 
as quantum mechanics, Yang-Mills type gauge theories, quantum non-linear optics, inte- 
grable systems and (quasi-)exactly solvable models, to name a few (see e.g. [1]-[10]). 

One of the reasons for their increasing prevalence stems from the realization that 
traditional linear Lie algebras describe only a very restrictive subset of linear symmetries 
and that many physical systems do in fact possess non-linear symmetries such as those in 
which the commutations of the symmetry algebra generators return polynomial terms. 

Due to their importance, a number of studies have been undertaken to investigate the 
mathematical properties of these algebras [11]. In particular, differential realizations of 



certain quadratic and cubic algebras have been explored in [12, 13] and also in [14, 15] in 
connection with the theory of quasi-exact integrability [16, 17, 18]. 

In this article, we introduce a novel class of higher order polynomial deformations 
of the classical A\ Lie algebra and construct their unitary representations in terms of 
boson operators and single-variable differential operators. We will then use the differential 
realizations of these algebras and the Functional Bethe Ansatz method (see e.g. [19, 20]) 
to obtain one of the main results of this paper, that is the exact eigenf unctions and energy 
eigenvalues of the following class of Hamiltonians 

2 2 

H = Y i w i N i + Y f w ij N i N j + g (a\ s a r 2 + a{a\ r ) , r, s e N, (1.1) 

where and throughout (aj) are boson or photon annihilation (creation) operators with 
frequencies Wi, iVj = ajaj are number operators, and g are real coupling constants. 
Without loss of generality, in the following we will identify w±2 with W2i- Hamiltonians 
(1.1) appear in the description of various physical systems of interest such as non-linear 
optics [8, 9] and Bose-Einstein Condensates (BECs) [21]- [24]. For instance, the non- 
diagonal terms in (1.1) describe processes of multi-photon scattering and higher-order 
harmonic generation in quantum nonlinear optics. Let us point out that (1.1) is a two- 
mode version of the more general multi-mode Hamiltonian considered in [25] in which the 
quasi-exact solvability of the multi-mode system was established by a different procedure 
and without giving exact solutions (see also [26] for the case of third-order harmonic 
generation)*. Hamiltonians for the special cases of s = r = 1 and s — 2, r = 1 have also 
been studied using Algebraic Bethe Ansatz (ABA) method [27]. 

This paper is organized as follows, in section 2 we propose a class of generalized 
polynomial su(l, 1) algebras and derive their boson realizations. In section 3, we use these 
deformed su(l, 1) algebras as base elements to generate higher order polynomial algebras 
via a Jordan-Schwinger like construction method. We then identify these algebras as 
the dynamical algebra of the Hamiltonian (1.1) in section 4 and solve for the eigenvalue 
problem in general via the Functional Bethe Ansatz method. In section 5, we present 
explicit results for the Hamiltonian (1.1) when r, s < 2 and r = s = 3. In section 6, we 
establish the spectral correspondence of these specific models with quasi-exactly solvable 
(QES) Schrodinger potentials. Finally we summarize our results in section 7 and discuss 
further avenues of investigation. 

*We became aware of these three references after submitting our work. We thank one referee for 
pointing them out. 
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2 Polynomial deformations of su(l,l) algebra 



Let k be a positive integer, k — 1, 2, • • •. We start off by proposing a class of polynomial 
algebras of degree k — 1 defined by the commutation relations 

[Qo,Q±] = ±Q±, 

[Q+,Q-\ = (fc) (Qo)-0 W (Qo-l), (2.1) 

where 

^(Q^-A(Q + i-J,) + fj(i^-I) (2.2) 

is a /c'^-order polynomial in Qo- The algebra admits Casimir operator of the following 
form 

c = q_q+ + 0«(q o ) = Q+Q- + w (Qo - 1). (2.3) 

For k — 1 and k = 2, (2.1) reduces to the oscillator and sw(l, 1) algebras, respectively. 
Thus, the algebra (2.1) can be viewed as polynomial extensions of the linear su(l, 1) and 
oscillator algebras. 

Similar to the su(l, 1) algebra case, unitary representations of (2.1) are infinite dimen- 
sional. In this section, we shall concentrate on the following one-mode boson realization 
of the algebra, 

Q+ =w (a,)t - Q - = m* (a)t ' Qo= H a,a+ ^)' (2 - 4) 

In this realization, the Casimir (2.3) takes the particular value, 

We now construct the unitary representations corresponding to the realization (2.4) 
in the Fock space Tib- There are k lowest weight states, 

|0>, ( fl t)|0> (at^lO). (2.6) 

Writing these lowest weight states as \q, 0) using the Bargmann index q, we have 

Qo\q,0) = q\q,0), Q_|g,0) = 0. (2.7) 

It follows from (2.3) and (2.5) that n?=i (<? + ~~ f) = ^' ^ rom wnicn we S et 

_ J_ fc + 1 2fc + 1 (fc- l)fc + 1 

^ ~~ £2 ' £2 ' £2 ' " ' ' £ 2 ' v2-o) 

This means that the boson realization (2.4) corresponds to the infinite dimensional unitary 

representation with particular q values (2.8). In other words, the %b decomposes into the 

i i (fc-i)fc+i 

direct sum T-L b = Hf 1 © • • • © H b ^ of k irreducible components %Y ■> ■■■■>%b ^ • 
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Noting kq—j: = 0, 1, .., k — 1 for all the allowed q values given in (2.8), we can write the 
lowest weight states (2.6) as \q, 0) = (a)) kq ~Ti |0). The general Fock states \q, n) ~ Q+\q, 0) 
in the irreducible representation space are then given by 



Ak(n+q- A ) 



\q,n) 



:|0). 



(n + g - p) 

It is easy to show that Q , Q± and C act on these states as follows 



Qo\q,n) 

Q+\q,n) 

Q-\q,n) 
C\q,n) 



= (q + n)\q,n), 
= Y[(n + q + 



ik-\" 



i=i 



n [ n +i 

fi-k 1 

i=i V 

0,1,- 



k 2 

i- l)k+l 



i=i 

k 

n 



A: 



A; 2 



A' 2 



\Q,n), 



q,n+ 1), 

|g,n- 1), 



1 k + 1 



9 = T?> 



A; 2 ' A; 2 ' 



(fc- l)fc + 1 
Ac 2 ' 



(2.9) 



(2.10) 



3 Polynomial algebras via Jordan-Schwinger type con- 
struction 

The unitary representations of the polynomial algebras discussed in the preceding section 
are all infinite dimensional. In this section, we shall employ a Jordan-Schwinger like 
construction [12, 13], to derive polynomial algebras that have finite dimensional unitary 
representations. Towards this end, we consider two mutually commuting polynomial 
algebras introduced in the preceding section, j<3+\ Q-\ <3o^} °f degree (ki — 1) and 
{Q ( +\ Q { -\Q { v ] } of degree (k 2 - 1), where k u k 2 = 1, 2, • • , 
Introduce new generators, 

Q+ = Q { +Q-\ Q-=Q?Q W , Qo = \{qS ) -QP). (3.1) 

We can easily show that Q ,± form a polynomial algebra of degree (ki + k 2 — 1) which 
close under the following commutation relations: 

[Qo,Q±] = ±Q±, 

[Q+,Q-] = ^ +k2 \Q ,C)-^ kl+k2 \Qo-i,C), (3.2) 

where 

£ = l(Q ( o ] +Q { o } ) (3-3) 
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is the central element of the algebra, 



[A Q±,o] = (3.4) 

and 

(So. C) = - ft (C + So + ± - i) ft (C - (& + 1) + 1 - 1) (3.5) 

is a (fci + A; 2 ) <?i -order polynomial in Qo and the central elements C. The Casimir operator 
of (3.2) is given by 

C = Q-Q + + ^ (kl+k2) {Qo,£) = Q+Q- + <p (fcl+fca) (Qo " 1, A- (3-6) 

For ki + k 2 = 2, i.e. k± — k 2 — 1, the polynomial algebra (3.2) reduces to the linear sw(2) 
algebra. So the algebras defined by (3.2) are polynomial deformations of su(2). 

In terms of two sets of mutually commuting boson operators acting on the tensor 
product of the Fock spaces, we have the realization (i — 1, 2 below) 

This realization gives rise to finite dimensional representations of the polynomial algebra 
(3.2). To show this, let \q±, ni), \q 2 , n 2 ) be the one-mode Fock states of the algebras 
{QoX} > {QoX} respectively, where m, n 2 = 0, 1, • • •, and q x = ^,^t 1 ,---, ( fcl ~^ fcl+1 
and g 2 = p-, ^pr~> • • • , ( fc2 ~^ fc2+1 . The representations of {<2o,±} corresponding to the 
realization (3.7) are then given by the two-mode Fock states \qi, ni)\q 2 , n 2 ). Since C 
is a central element of the algebra, it must be a constant, denoted as I below, on any 
irreducible representations. This imposes a constraint on the values of n\ and n 2 , 

£\qi,ni)\q 2 ,n 2 ) = ^(q 1 + m + q 2 + n 2 )\qi, rii) \q 2 , n 2 ) 

= i|gi,ni)|g2,w 2 ). (3.8) 

That is rii + n 2 = 21 — (qi + q 2 ). Thus obviously 21 — qi — q 2 take only positive integer 
values, i.e. 

2Z - gi - g 2 = 0, 1, • • • ■ (3-9) 
It follows that the Fock states corresponding to the realization (3.7) are 

\qi,q 2 ,n,l) = \qi, n)\q 2 , 21 - q x - q 2 - n) 



t kAn+ qi -X\ . k 2 hl-qi-n-X 

(a[) v k i J (a 2 ) v 2 



vi ~ n 



= |0>, 



(k 1 (n + q 1 -^))\ ] /(k 2 (2l-q 1 -^ 
n = 0,l,---,2l- qi -q 2 , (3.10) 
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noting that 21 — qi — q 2 is always less than or equal to 21 — q-i — ^. This gives us the 
21 — qi — q 2 + 1 dimensional irreducible representation of (3.2), 

Qoki,92,w,Z) = (gi - Z + n)|gi,g 2 ,w,Z), 



Q+|gi,g 2 ,n,Z) = JJ 2Z - qx 



n 



k 2 {i-i) + r 



h 



2 



kl ( jki — 1 \ 1 

xJJ[n + gi + — 72 — \qi,q2,n + l,l), 



Q-\qi,q2,n,l) = JJ ( 2Z - q ± - n + 



j=l \ ft l / 

k2 f iA: 2 — 1\ * 

i \ k 2 J 

x ft ( n + 9i - (j ~l ) 2 A:i + 1 l 2 92, n - 1, Z). (3.11) 



By using the Fock-Bargmann correspondence, 

_d_ 



d z ni 
a\ — > Zi, Oi — > — , K) — >■ -^==, (3.12) 



we can make the following association 

fcitn+gi-i) fe 2 (2i-gi-n-i) 

1 2 

/2 z 

\qi,q 2 ,n,l) — » -= 1 =-^= (3.13) 

'(Mn + ft-^y^Z-ft-^-n))! 

Now since l,qi,q 2 ,k±,k 2 are constants, we can map the states |?i,?2,^, above to the 
monomials in z = Zi 1 /- 2 ^ 2 , 



^qi,(j2,n,/(^) 



(h(n + q ± - ^))\Mk 2 (2l - qi -^- n ))\ 



n = 0,l,---,2l- qi -q 2 . (3.14) 

The corresponding single-variable differential operator realization of (3.2) takes the fol- 
lowing form 

d 

Qo = z— + q 1 - Z, 

dz 

These differential operators form the same 21 — qi — q 2 + 1 dimensional representations in 
the space of polynomials as those realized by (3.7) in the corresponding Fock space. We 
remark that because HjLi (<?i — ^~ 1 fc ^ 1+1 ) = for all the allowed qi values there is no z^ 1 
term in Q_ above and thus the differential operator expressions (3.15) are non-singular. 
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4 Exact solution of the 2-mode boson systems 

We now use the differential operator realization (3.15) to exactly solve the 2-mode boson 
Hamiltonian (1.1). 

By means of the Jordan- Schwinger type construction (3.1) and the realization (3.7), 
identifying h\ with s and h 2 with r, we may express the Hamiltonian (1.1) in terms of the 
generators of the polynomial algebra (3.2), 

2 2 

H = Y^ w i N i + J2 w a N i N 3 + g^^(Q+ + Q-) (4.1) 

with the number operators having the following expressions in Q and C 

N l = s(Q + C)--, N 2 =r(C-Q )--. (4.2) 

S T' 

Keep in mind that {Q Ti o} in (4.1) as realized by (3.7) (and (3.1)) form the (2/ — q\ — 92) + 1 
dimensional representation of the polynomial algebra (3.2). This representation is also 
realized by the differential operators (3.15) acting on the (21 — q 1 — q 2 ) + 1 dimensional 
space of polynomials with basis jl, z, z 2 , z 2l ~ qi ~ q2 }. We can thus equivalently represent 
(4.1) (i.e. (1.1)) as the single-variable differential operator of order max{s,r, 2}, 

2 2 

H = J2 W i N i + Y. W V N i N j 

i i,j 

iV ' = S <4 + "»-? JV> = r<2l -i. (4.4) 
We will now solve for the Hamiltonian equation 

Hijj(z) = EiP(z) (4.5) 

by using the Functional Bethe Ansatz method, where ip(z) is the eigenfunction and E is 
the corresponding eigenvalue. It is easy to verify 

Hz m = z m+1 gj\r [2I — q\ — — ^— rnj + lower order terms, m G Z + . (4.6) 

This means that the differential operator (4.3) is not exactly solvable. However, it is quasi 
exactly solvable, since it has an invariant polynomial subspace of degree (2/ — q 1 — q 2 ) + 1: 

HVCV, V = span{l,z,...,z 2 '- 91-92 }, dimV = 21 - qi - q 2 + 1. (4.7) 



with 
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This is easily seen from the fact that when m = 21 — q 1 — q 2 the first term on the r.h.s. 
of (4.6) becomes z 2l ~ qi ~ q2+1 g JYj=i r (?2 — ^~^ r+1 ) which vanishes identically for all the 
allowed q 2 values. We remark that the quasi-exact solvability of the system is connected 
with its quantum integrability, i.e. with the fact that there exists quantum operator 
coinciding with a linear combination of the operators Aq and N 2 which commutes with 
the Hamiltonian (1.1). 

As (4.3) is a quasi exactly solvable differential operator preserving V, up to an overall 
factor, its eigenfunctions have the form, 



M 



i>(z) = J] (z - at) , (4.8) 



i=i 



where M = 21 — q± — q 2 (= 0, 1, • • •), and {ctj \ % — 1, 2, • • • , M} are roots of the polynomial 
which will be specified later by the associated Bethe ansatz equations (4.14) below. We 
can rewrite the Hamiltonian (4.3) as 

max{r,S,2} / 7 \ i 

H= E P ^[d- Z ) + P oW (4-9) 



1=1 



where 



P (z) = zgflr(2l- qi - ^ 1 j r + 1 > 



i=i 



+wu (sqi ~ ^) + w 22 (r (2/ - qi) - ^ 
+2w 12 (s qi - -) (r (21 - Ql ) - - 



+w 1 ^sq 1 - ] -^+w 2 (r(2l-q 1 )-^j (4.10) 
and Pi{z) are the coefficients in front of d l /dz l in the expansion of (4.3) (see the Appendix), 

Pi(z) = 9s s z^±( ± f[ WUjm 

k=i \h<...<l k j^h^-^l k J 

+g(- r yz^±( ± n b\l Ki 

k=i \h<...<l k j^h^-^l k / 

+F5 i:2 z 2 + D6i, lZ . (4.11) 



In the above expression, 



A = qi 



[i-l)s + l 
s 2 



Lk,k — 1, 
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k-l 



L k ,i = E n i( n 2 - l)--(^fe-i — {k — i) + 1), i < k, 

ni<...<n fe _i 

F = w 22 r 2 + wns 2 - 2wi 2 sr, 

D = w 22 r 2 (l - 2(2/ - gi - — )) + w n s 2 (l + 2( gi - -) 

+2u?i 2 rs ^2(Z - 91) + - - 1^ + wis - w 2 r. 
Dividing the Hamiltonian equation Hip = Eip over by tp gives us 

1 

/,. h<...- h (z-a h )...{z-a h 



E 



TT I max{r,S,2} M 

= E E 



(4.12) 



(4.13) 



The l.h.s. of (4.13) is a constant, while the r.h.s is a meromorphic function in z with at 
most simple poles. For them to be equal, we need to eliminate all singularities on the 
r.h.s of (4.13). We may achieve this by demanding that the residues of the simple poles, 
z — oti, i — 1, 2, M should all vanish. This leads to the Bethe ansatz equations for the 
roots {a.j} : 



,x{r,s,2} 

E 



M 

E 



Pi(a p )i\ 



+ P 1 (a p ) = 0, p=l,2,--.,M. (4.14) 



i=2 h<l 2 <-<li-i^p ( a P a h)'"\ a P a h-U 

The wavefunction ip(z) (4.8) becomes the eigenfunction of H (4.3) in the space V provided 
that the roots {«;} of the polynomial ip(z) (4.8) are the solutions of (4.14). 

Some remarks are in order. It is easily seen (from (4.8) and (4.15) below) that Hip/ip is 
regular at z — ±oo. When (4.14) is satisfied, the r.h.s. of (4.13) is analytic everywhere in 
the whole complex plane and thus must be a constant by the Liouville theorem. Therefore 
the Bethe ansatz equation (4.14) is not only necessary but also sufficient condition for the 
r.h.s. of (4.13) to be independent of z. 

To get the corresponding eigenvalue E, we consider the leading order expansion of 



M 



„M „Af-l 

Z — Z 



+ 



i=l 



It is easy to show that Q±fli J ( z ) have the expansions, 



Q-ip = z 



—z 



M 



nk+i 

j=i \ 



M-liVf)_ 



n 2/ 



Q2 



(j-l)r + r 



(j-l)s + l 



M 

E«i 
i=i 



+ 



(4.15) 



Qoifj = z M (l-q 2 ) + ---. 

Substituting these expressions into the Hamiltonian equation (4.5) and equating the z M 
terms, we arrive at 



E = w u (s(2l - q 2 ) - ^ + w 22 (rq 2 



;) 2 
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+2w 12 (s(2l - q 2 ) - (rq 2 ~ J 
+wi (s(2l - q 2 ) - ^ + w 2 (rq 2 - ^ 



~9 



U r \Q2 + l ~ 2 



J'=l 



M 



E«i. ( 4 - 16 ) 



where {aj} satisfy the Bethe ansatz equations (4.14). This gives the eigenvalue of the 
2-mode boson Hamiltonian (1.1) with the corresponding eigenfunction ip(z) (4.8). 

5 Explicit examples corresponding to BECs 

We will now work out in complete detail the Bethe ansatz equations and energy eigenval- 
ues of the Hamiltonian (1.1) for the special cases of s, r < 2 and r = s = 3. These models 
arise in the description of Josephson tunneling effects and atom-molecule conversion pro- 
cesses in the context of BECs. 

A. s = 1, r = 1 

The Hamiltonian is 

2 2 

H = J2 WiNi + WijNiNj + g (a\a 2 + ai a\) . (5.1) 

This is the so-called two coupled BEC model and has been solved in [27] via a different 
method, i.e. the ABA method. From the general results in the preceding section, in this 
case, we have q\ = q 2 = 1, which means that 21 — q x — q 2 = 2(1 — 1) = 0, 1, • • •. That is 
I — 1 = 0, |, 1, • • •. The differential operator representation of the Hamiltonian (5.1) is 

H = P*{z)^ + Pi{z)± + P Q {z), (5.2) 



where 



with 



P 2 (z) = A n z 2 , 

P 1 (z) = -gz 2 + B n z + g, 

P (z) = 2(l-l)gz + D n (5.3) 



Au = w u + w 22 - 2w l2 ^ 0, 

B n = w l -w 2 + wu + (5 - 41) w 22 + (41 - 6) w 12 , 

D n = 2(1 - l)w 2 + 4(1 - l) 2 w 22 . (5.4) 
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The Bethe ansatz equations are given by 

2( ^ 1} 2 g + gnap - ga 2 p 
^ «i-« P Ana* 

and the energy eigenvalues are 

2(1-1) 

E = 4w u (l ~ I) 2 + 2 Wl (l -l)-g E «i- ( 5 - 6 ) 

i=i 

B. s = 2,r = 1 

The Hamiltonian is 

2 2 

H = ]T lUiiVi + £ WijNiNj + g (a ]2 a 2 + a?4) . (5.7) 

i i,3 

This is the homo-atomic-molecular BEC model and has been solved by the ABA method 
[27]. Specializing the general results in the preceding section to this case, we have qi = \ 
or |, and q 2 — 1. The differential operator representation of the Hamiltonian (5.7) is thus 

H = P 2 (z)-^ + P 1 (z)± + P (z) (5.8) 



where 



P 2 (z) = A 21 z 2 + 4gz, 

P 1 (z) = -gz 2 + B 21 z + 8gq 1 , 

P (z) = g (2l- qi -l)z + D 21 (5.9) 



with 



A 21 = Awn + w 22 - Aw 12 , 

B 21 = 2w 1 -w 2 + 2w 11 (l + Aq 1 )+w 22 (3 + 2q 1 ~Al) 

+w 12 (-7 -8q! + 8l) , 
D 21 = 2w 1 (q 1 --^+w 2 (2l-q 1 -l) 

+Aw u ~\) + w 2 2 (2/ - 1 - qif 

+ Aw l2 [q 1 -^j{2l-l-q 1 ). (5.10) 

The Bethe ansatz equations are 

2l X- qi 2 8gq 1 + B 21 a p -gal 

1^ = U rT\ i p = l,2,---,2Z-l-gi 5.11 

r£ ai-a p a p (A 21 a p + Ag) 
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and the energy eigenvalues are given by 

E = 2w 1 {2l--J+Aw 11 (2l--J -9 a i- ( 5 - 12 ) 



C. s = 2,r = 2 

The Hamiltonian is 

2 2 

if = 2 ^ + E tfij^i^ + 5 (a\ 2 al + afc?) . (5. 13) 

i i,j 

This gives another model of the atom-molecule BECs. To our knowledge, this model has 
not been exactly solved previously. Applying the general results in the preceding section, 
we have in this case qi — \, § and q 2 — \, §■ The differential operator representation of 
the Hamiltonian (5.13) is 

H = P 2 (z)-^ + P 1 (z)^ + P (z), (5.14) 

where 

P 2 (z) = Agz 3 + AA 22 z 2 + Agz, 
Pi(z) = B 22 z 2 + D 22 z + Sgq 1: 

p ( z ) = F 22 z + G 22 , (5.15) 

with 

-422 = wii + w 22 - 2wi 2 , 
B 22 = 8^(l + gi -2Z), 

D 22 = 2w 1 -2w 2 + 2wn(l + 4g 1 ) + 2w 2 2(3-8/ + 4g 1 ) 

+Sw 12 (-1 -2gi + 2Z), 
F 22 = Ag(2l- qi --^ (2l- qi -^j, 

G 22 = 2 Wl ( qi -\)+ 2w 2 ( 2 1 - qi - \) + Aw u ( qi - \) 

2 



4/ V 4/ V 4 



+8w;i2 (ft - ^) (2/ - ft - ^ + 4w 22 (2/ - ?i - J) • (5-16) 

Note that (21 - ft - q 2 )(2l - ft + q 2 - 1) = (2Z - ft - 1/4)(2Z - ft - 3/4) for g 2 = 1/4,3/4. 
The Bethe ansatz equations read 

= 1 — 7 — 2— 1 r^' P=l,2,---,2Z-gi-g 2 (5.17) 
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and the energy eigenvalues are 



E = Aw u ( 21 - q 2 - ^ + Aw 22 (ft - | 



+8u> 12 (21 - q 2 - ^ (ft - J) + (2/ - ?2 - J) 

/ 1\ / 1\/ 3\ 2l ~ qi ~ q2 

+2w 2 [q 2 --j-Ag[q 2 + -j(q 2 + -j £ a,. (5.18) 



D. s = 3, r 



The considered examples with s, r < 2 may in principle be treated using the ABA 
method based on Lie algebra su(2) (without any polynomial deformations). We now 
present an explicit example for which the ABA method is not applicable . The Hamilto- 
nian is 

2 2 

H = WiNi + WijN i N j + g (afa\ + afaf.) . (5. 19) 

This is a non-linear optical model with third-order harmonic generation. Specializing the 
general results in the preceding section to this case, we have ft, ft = §, §, or |. The 
differential operator representation of the Hamiltonian (5.19) is 

, N d 3 , N d 2 , N d / \ , N 

ff = P3W ^ + P2W ^ + Pl{z) Tz + PoW ' (5 - 20) 

where 

P 3 (z) = 27g(-z 4 + z 2 ) 

P2(z) = A 33 z 3 + B 33 z 2 + D 33 z, 

P 1 {z) = F 33 z 2 + G 33 z + K 33 

P {z) = R 33 z + S 33 (5.21) 

with 

A 33 = 9g (18/ -9ft -13), 
5,33 = 9(w n + w 22 - 2w 12 ), 
D 33 = 9<?(9ft + 5), 

F 33 = 9g (-36/ 2 - — + 34/ + 36/ft - 9ft 2 - 17ft J , 

G33 = 3wi - 3u> 2 + wn(7+ 18ft) + 2wi 2 (-9- 18ft + 18/) + w 22 (7 + 18ft), 
^33 = ^(qi + Qql + l), 

i?33 = 27,(2/-ft-i)( 2 /-ft-l)( 2 /-ft-I), 

S 33 = 9W U (ft - ^) + 9W 22 (2/ - ft - + I8W12 (ft - -) (2/ - ft - 

+3 Wl (ft - g) + 3^ 2 (2/ - ft - ^ . (5.22) 
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The Bethe ansatz equations read 

2l ~^ 162^ - a 2) + 2( ^ 33a 3 + ^ 33a 2 + D33Qp) 

= F 33 a 2 p + G 33 a p + K 33 , p = 1, 2, • • • , 21 - q ± - q 2 (5.23) 
and the energy eigenvalues are 

E = 9w n (21 -cb-^J + 9w 2 2 (<?2 - 



+ 18w 12 (21 - 52 - ^) (?2 3u, i ( 2/ - <?2 - £ 

+3 W 2(g2-gj-27^g 2 + -j(g 2 + -j(g 2 + -j g a, (5.24) 



6 Spectral Equivalence with QES Schrodinger Poten- 
tials 

The Hamiltonians in section 5 correspond to second order differential operators and can 
be mapped to Schrodinger equations with QES potentials via a suitable similarity trans- 
formation and change of variables [28]. 

Explicitly, if H is written in the following form 

H = P{z)^ + (q(z) + Ip'(z)) 4- + R(z) (6.1) 



dz 2 V w 2 v 7 dz 
then it can be mapped to a Schrodinger operator, 

H = -e w ^He~ w ^ = ~ + V(x), (6.2) 

dx z 

where the variable x and z are related by (we assume z — z(x) is invertible on a certain 
interval to give x = x(z) ) [29], 

x = x(z) = ± f -fi= (6.3) 

and W(x) is given as 

The potential function is given by 



(6.5) 

z=z(x) 
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Then the solutions of the 2nd order ODE Hip(z) = Etp(z) with eigenvalue E is mapped 
to solutions of the Schrodinger equation 



Hip(x) = Eip(x) (6.6) 
with eigenvalue E = —E and corresponding Schrodinger wavefunction 

$(x) =e- w{x) i>{z(x)). (6.7) 

We will not discuss the square integrability of the Schrodinger wavefunction if)(x), 
but derive the explicit Schrodinger potentials corresponding to the special models in the 
preceding section. 

I. s = l,r = 1: 

For this case we have 

P(z) = A u z 2 , 

Q(z) = -gz 2 + (B 11 -A 11 )z + g, 

R(z) = 2(1 - l)gz + D u . (6.8) 

From (6.3), we obtain 

z {x) = e^ x (6.9) 

The potential is 

T/ , v {gz 2 + {A u - B n )z - g) (gz 2 + {3A U - B u )z - g) 
V[X) ~ AA llZ * 

-g(2l-l)z + ^^-D n 
9 2 , f n B 



2A 



cosh(2 Ja^ix) +g(2- sinh(JA 1 ~ 1 x) 
i 1 V An / 1 



-(2/ - l)ge^* - D u + ( A n ~ Bn) 2 ~ 2g\ (g 1Q) 

II. s = 2,r = 1: 

In this case we have, 

P(z) = A 2 iz 2 + 4gz, 

Q(z) = -gz 2 + (B 21 - A 2 i)z + 2g(Aqi - 1), 

R(z) = g(2l-qi-l)z + D 2 i (6.11) 

and qi = \ or |. From (6.3), we derive 

z(x) = ^- (cosh( y/A^x) - 1) . (6.12) 
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The potential is 



V(x) 



(gz 2 + (An - B 2l )z + 2g - gggQ (ffz 2 + (3A 21 - B 21 )z + 6ff - 8ggi) 

4(A 21 z* + Agz) 



(01 \ i ^21 — A 21 
-g[2l - qi)z H D 21 



11 

A? 



tanh 2 ( — ' ■ 



,r | siiilr I — x 



21 



V . , 2 / V^l" 

— — sinn 



x ) + 2A 2 i - B 21 



121 



4A 



21 



(A 21 - B 21 ){3A 21 - B 2l ) + 8(? 2 (1 - 2gi) 



tanh* 



A 



X 



21 



, (3 - 8 qi )A 21 + (4gi - 2)^21 5 2 i - A 2 i n 

H o / AT— \ 1 n -^21- 



4 cosh 2 (^x) 

Here we have used (1 — 4gi)(3 — Aq{) = for the two allowed q± values q\ = \ or |. 
Let us consider the special case of A 2 i = 0. In this case, 



(6.13) 



z(x) = gx 2 



(6.14) 



as can be seen from the A 21 — > limit of (6.12). The Schrodinger potential (6.13) then 
becomes 

Vt \ ^ 6 o 4 ^ ^+8^(1-4/) 2 . p n 

V^(x) = — :r D - — B 21 x H — x +qiB 21 -D 21 . 

16 8 16 

This is a non-singular sextic potential. 



(6.15) 



III. s = 2,r = 2: 

For this case we have 



P(z) = Agz 3 + AA 22 z 2 + Agz, 

Q(z) = B 22 z 2 + {D 22 - AA 22 )z + 8g qi - 2g, 

R(z) = F 22 z + G 22 



(6.16) 



and qi = \ or |. From (6.3), we obtain 



z( x ) = g 3 p(gnx;g 2 ,g 3 ) 



A 22 

3<? : 



(6.17) 



where p(x;g 2 ,g 3 ) is Weierstrass's elliptic function with invariants g 2 and g 3 given by 



4 ^ (AL \ A ( 2A 2 \ 

3^(^-3), ^=27^22(9-^). (6.18) 
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Hereafter we will denote p(x; g 2 , g%) simply as p(x). The potential is computed as follows 

V(x) = (b 22 z 2 + (D 22 - AA 22 )z + 8g qi - 2g) 

(B 22 - \2g)z 2 + (D 22 - \2A 22 )z + 8g qi - 6g 



16 (gz 3 + A 22 z 2 + gz) 

D 22 - AA 22 
+ (B 22 — F 22 )z H G 22 

= E c *~ A „ 1 s , +(B 22 -F 22 )g sp(g*x) 



i=i 



+ ^ + 2 G22 ' (6 - 19) 



where 



ci = 2( ? (4g 1 -3)( J D 2 2-4A 2 2) + 2^(4^-1) (D22-I2A22), 

c 2 = 2g(A qi - 3)B 22 + 2g{A qi - 1){B 22 - I2g) + {D 22 - AA 22 ){D 22 - 12A 22 ), 

c 3 = B 22 (D 22 -12A 22 ) + (B 22 -12g)(D 22 -AA 22 ), 

c 4 = B 22 (B 22 - 12g). (6.20) 

Here we have used 4(gz 3 + A 22 z 2 + gz) = Ap(g^x) 3 — g 2 p(g^x) — g^ = p'(g^x) 2 and 
(1 — 4gi)(3 — 4^) = for the two allowed qi values q ± — \ or |. 



7 Discussion 

Let us now quickly summarize the work. We began by constructing the boson represen- 
tation of a class of su(l, 1) polynomially deformed algebras (2.1), deriving their infinite 
dimensional Fock space realization and lowest weight state parametrization. We then 
used the Jordan-Schwinger like construction to get the polynomial algebra (3.2) which 
possesses finite dimensional irreducible representations. We used the differential realiza- 
tion of (3.2) to rewrite the Hamiltonian (1.1) as QES differential operators acting on 
the finite dimensional monomial space. The exact eigenf unctions and eigenvalues of the 
Hamiltonian were then found by employing the Functional Bethe Ansatz technique. As 
examples, we provided some explicit expressions for the BEC models which correspond 
to the r, s < 2 cases of (1.1) and established the spectral correspondence of these specific 
models with QES Schrodinger potentials. 

In deriving our results, we showed that in general the Hamiltonians defined in (1.1) 
are QES differential operator of order 3 or higher. This paper provides an algebraization 
of such higher order QES differential operators and unravels the dynamical polynomial 
algebra symmetry of (1.1). It also shows that the Functional Bethe Ansatz method 
provides a simple way to find exact eigenvalues and eigenfunctions of such higher order 
differential operators. 



17 



There are a number of extensions that we plan to pursue in this line of investigation. 
Firstly, we note that the Jordan- Schwinger like construction can be extended straightfor- 
wardly to study other non-linear quantum optical models such as the general multi-mode 
boson Hamiltonians of the form 

k+k' k+k' 

H = Y. "•'• V < • E w ij N i N j 

+ 9 (a\^ ■ ■ ■ 4 mfe cr • • • cr + <c • • • «r4?i +i • • • 4:r ) • (7.1) 

Results on on this and other models of physical interest will be presented elsewhere. 

Acknowledgments: This work was supported by the Australian Research Council. The 
authors would like to thank Ryu Sasaki for very valuable comments and suggestions which 
lead to significant improvement of the presentation of the paper. 

8 Appendix 

In this appendix, we work out the expansion coefficients in front of ^ in the xpansion of 

n£i0±+4). 

First, we see 



nfe+^ = n*+ eiia \z-+ e n a ^ 

i=l \ az / i=l Vi=li#ii / aZ V/r / \ a2 7 



+•••+( e n a 



7 \ "I 

d \ 



Z dz 



(8.1) 



and since 



/ d \ k d k ( k ~ x \ d k ~ l ( fc_1 \ d k ~ 2 

1=1 

where 

^it,fe = 1 

fe-i 

L fcji = n i( n 2 - l)...(n fc _j - (fc - i) + 1), i < k (8.3) 

ni<...<n fc _i 

We can regroup the equation (8.1) as 

8(*= + *) 



18 



m I I m m \ / m m \ \ i 

=n^+k En^ +-+w e n ^ 4 

»=1 V \il = l^il / yi<J2<-"<jmi#l"-#m // 

^ 2 e n a e n ^ 4 

\h<32i^h^32 J \3l<32<-<3mi^3v¥3m ) ) 

+ higher order terms 

= n^+EE e n A A L ^ l h- (8 - 4) 
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